A solvable quantum antiferromagnet model 

Bikas K. Chakrabarti If and Jun-ichi Inoue f 
I Saha Institute of Nuclear Physics, 1/AF Bidhannagar, Kolkata 700064, India. 

bikask. chakrabartiOsaha. ac . in 
t Graduate School of Information Science and Technology, Hokkaido University, 
N14-W9, Kita-ku, Sapporo 060-0814, Japan. 
j_inoue@complex . eng .hokudai .ac.jp 

Abstract 

We introduce a quantum antiferromagnet model, having exactly soluble thermodynamic proper- 
ties. It is an infinite range antiferromagnetic Ising model put in a transverse field. The free energy 
gives the ground state energy in the zero temperature limit and it also gives the low temperature 
behaviour of the specific heat, the exponential variation of which gives the precise gap magnitude in 
the excitation spectrum of the system. The detailed behaviour of the (random sublattice) staggered 
magnetisation and susceptibilities are obtained and studied near the Neel temperature and the zero 
temperature quantum critical point. 



With the realisation in the mid last century, that the Neel state can not be the ground state, or even an 
eigenstate, of the quantum Heisenberg antiferromagnets, considerable efforts have gone in search for the 
nature of the ground state for such quantum antiferromagnets. Searches have also been made for low lying 
eigen states of Heisenberg antiferromagnets even in the semi-classical limits, or in some other variants 
of the model where exact ground state and other eigenstates are available [1,2]. In particular, the exact 
dimerised (two-fold degenerate; disordered) ground state in the one dimensional next nearest neighbour 
interacting Heisenberg antiferromagnet (having strength ratio 1/2 between the next nearest and nearest 
neighbours) has been found out [3] and this observation has attracted considerable attention in the 
context of high temperature superconductors occurring in materials having antiferromagnetic properties 
[1,2]. Regarding the excited states in such quantum antiferromagnets, the Haldane conjecture [2,4] states 
that the integer spin systems are massive (have a gap), while the half-integer ones are massless. 

The low temperature limit of the exact free energy obtained here for the (spin-i) long-range interacting 
transverse Ising model shows that the specific heat has an exponential (in gap and inverse temperature) 
variation, giving the precise magnitude of the gap in the model. The order-disorder transition, driven 
both by temperature and the tunneling or transverse field, are investigated studying the ordering and 
susceptibility behaviours. 

The model we study here has the Hamiltonian 

H^j^J2S'S^~hJ2St-Tj2s:, (1) 

i.j i i 

where J denotes the long-range antiferromagnetic (J > 0) exchange constant and 5^ and denote the 
X and z components of the N Pauli spins {S = 1/2): 

^f - ( J J )' - ( J )' * = 1,2,. ..,7V. 

Denoting half of the randomly chosen lattice sites as members of sublattice A and the rest of B, and 
expressing the cooperative term in Hamiltonian (1) as the difference of two quadratics, the Hamiltonian 
can be rewritten as 
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In fact, the Hamiltonian (1) is meaningful when the spins in the cooperative term there belongs to 
two different (though randomly defined) sublattices, each consisting of N/2 spins. Using the Hubbard- 
Stratonovich transformation, the partition function can be expressed as 



Z = Tr 



Sa, Se 
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where f3 denotes the inverse temperature T and m+ 
magnetisation M 



are related to the (uniform; ^-component) 



^ and Mg of the sublattices A and B respectively: 



Ml,) and 



(M\ — Mg). The trace in (2) can now be easily performed [5] as each spin S can now be imagined to be 
isolated and present only in a vector field h having both z and x components: — T and /i^ ^ f3Jm. 
In the N oo limit, the free energy / per spin can be obtained from 
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where the sublattice magnetisations M\ and M^ are given by the (self-consistent) saddle point equations: 



Ml = 



Mi 



(-2 JM| + /i) 
v/(-2 JM| + /i)2 + r2 

(-2 JMl + /t) 
^(-2 JMl + /i)2 + r2 



tanh /3 J (-2 JM| + /i)2 + r2 



tanh/3./ (-2 JMl + h)^ + T^. 



(4) 



The spontaneous sublattice order Ml or Mg vanishes at the Neel phase boundary Tj\[{T); see Fig-Q] 
It may be noted that the choice of sublattices here is purely random and the sublattice order therefore 
looks like that of a glass. Also, deep inside the antiferromagnetic phase (at /3 ~> oOjT 0,h — 0), 
Ml = 1 = —Mg, so that the free energy / can be expressed as / ~ (1//3) log[l + exp(— 2/3A)] and the 
specific heat d^f/dT^ will have a variation like exp[— 2/3A(r)], similar to those of a two level system with 
a gap A(r) = v4j2 + r2. This is the exact magnitude of the gap in the magnon spectrum of this long 
range transverse Ising antiferromagnet. Actually, the above result for the gap can be seen directly from 

a b 



the effective Hamiltonian in (2) having a simple 



b -c 



matrix form where a 



2l3Jm and b = T, 
r, where the magnetisation 



giving the eigen values ±(3^/ A,Pm? + r2. This gives the gap /3A at T = 
m becomes equal to unity. 

In the classical case, near the Neel temperature = r/v(r = 0), we expand the equations (4) around 



Ma ~ 0, Mb ~ for ft. = and obtain M% 



-2J(3Mli + 0{{MIj)^),MI 



The only possible solutions for these linearized equations are Ml 



-M? 



-2J/3M1 + 0((M1)3). 
M^ = and Ml = 



-M| = M^ 



(T — T^)^/2^ where = 2 J. The longitudinal linear susceptibilities Xa/b are given by 



XA/B = lim/j^o dM%ig/dh = /3(1 - 2,/xBM)cosh /3(~2 JM^/^). Hence, x = Xa + Xb behaves as 



X 



T% + Tcos\y{2JM^)/T' 



(5) 



This gives x = 2/(r + T^) for high temperatures T > T°, (where Ma = -Mg = M'^ = 0) and at 
it grows upto a value 1/2 J and drops down again to as T — > 0, giving the well known cusp behaviour 
(see e.g., [6]) for the classical antiferromagnets. 

In the pure quantum case (T = 0, F 7^ 0), the tanh term in (4) equals to unity and with a similar 
expansion near the quantum critical point F^r = 2 J, one gets the longitudinal susceptibility x — Xa + Xb] 
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Figure 1: Variations of longitudinal sublattice magnetization and susceptibilities in the classical and quantum 
cases {J = 1). 



XA/B = dMX/g/dh ^T^il- 2 Jxbm)/[(2 JM^/^)^ + r^]^/\ We then have 

^ " [(rjvM^)2 + r2]3/2 + rjvr2' 

giving X — 2/(r + rjv) for T > Tn and growing upto a value 1/2 J at F = Fat and then decreasing 
eventually to a value x — F2/4J'^ as F — > 0. Again a cusp behaviour is seen for the susceptibility at 
the quantum critical point Fjv for such a quantum antiferromagnet (see Fig. The behaviour is of 
course qualitatively similar to that the classical critical (Neel) point. One can also study the transverse 
susceptibility x^{— dM^/dT) behaviour. One finds X;4 = = for F > Fjy and Xa ~ 1/2 J = Xb 
the ordered phase (see Fig.|21). 

In view of the intriguing ground state properties of quantum antiferromagnets [2] , indications coming 
often only from approximate theories (see e.g., [7] for a long range quantum Heisenberg antiferromagnet) 
or numerical simulations (see e.g., [8] and references therein), our proposed quantum antiferromagnetic 
model, having exactly soluble thermodynamic properties, should be of some interest. The model consists 
of an infinite range antiferromagnetic Ising system, put in a transverse field. The classical ground state 
of the model is highly degenerate. Although no signature of slow dynamics, like in glasses, can be seen 
here, the ordered state in the system corresponds to (quantum) glass-like system as well. The number 
of degenerate states can be estimated to be 0{2^^'^), which is larger than that for the Sherrington- 
Kirkpatrick model (~ 0('2O-28743JV'j^ jgj^ This may be compared and contrasted with the transverse Ising 
antiferromagnets on topologically frustrated triangular lattices studied extensively in the last few years 
[10] . The free energy in (3) gives the ground state energy in the zero temperature limit and it also gives the 
low temperature behaviour of the specific heat, the exponential variation of which gives the precise gap 
magnitude A(= V 4J2 + F^) in the excitation spectrum of the system. It may be noted that although it is 
a spin-1/2 system, because of the restricted (Ising) symmetry and the infinite dimensionality (long range 
interaction) involved, there need not be any conflict with the Haldane conjecture. Although our entire 
analysis has been for spin-1/2 (Ising) case, because of the reduction of the effective Hamiltonian in (2) to 
that of a single spin in an effective vector field, the results can be easily generalised for higher values of the 
spin S. No qualitative change is observed. The order-disorder transition in the model can be driven both 
by thermal fiuctuations (increasing T) or by the quantum fiuctuations (increasing F). These transitions in 
the model have been investigated here studying the behaviours of the (random sublattice) magnetisation 
and the (longitudinal and transverse) susceptibilities. No quantum phase transition, where the gap A 
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Figure 2: Variations of transverse magnetisation and susceptibility in the quantum case (J = 1). 



vanishes, is observed in the model, unhke in the one dimensional transverse Ising antiferromagnets [1,5]. 



Acknowledgements: The work was financially supported by IKETANI Science & Technology Foun- 
dation grant no. 0174004-D. One of the authors (B. K. C.) thanks University of Tokyo and Hokkaido 
University for their hospitality. We are also grateful to S. M. Bhattacharjee, I. Bose and Y. Sudhakar for 
illuminating discussions. 



References 

[1] See e.g., S. Sachdev, Quantum Phase Transitions, Cambridge Univ. Press, Cambridge (1999). 

[2] Sec e.g., E. Fradkin, Field Theories of Condensed Matter Systems, Addison- Wesley, Redwood (1991); 
I. Bose, in Frontiers in Condensed Matter Physics: 75th Year Special Publication of Indian Journal 
of Physics, Vol. 5, Eds. J. K. Bhattacharjee and B. K. Chakrabarti, Allied Publishers, New Delhi 
(2005). 

[3] C. K. Majumdar and D. K. Ghosh, J. Math. Phys. 10 1388, 1399 (1969). 

[4] F. D. M. Haldanc, Phys. Lett. A 93 464 (1983); Phys. Rev. Lett. 50 1153 (1983); I. Affleck, T. 
Kennedy, E. H. Lieb and H. Tasaki, Phys. Rev. Lett. 59 799 (1987). 

[5] B. K. Chakrabarti, A. Dutta and P. Sen, Quantum Ising Phases and Transitions in Transverse Ising 

Models, Springer, Heidelberg (1996). 

[6] C. Kittel, Introduction to Solid State Physics, John Wiley & Sons Inc., N. Y. (1966). 

[7] C. Kaiser and L Peschel, J. Phys. A : Math. Gen. 22 4257 (1989). 

[8] T. Roscilde, P. Verrucchi, A. Fubini, S. Haas and V. Tognetti, Phys. Rev. Lett. 94, 147208 (2005). 

[9] F. Tanaka and S. F. Edwards, J. Phys. F: Metal Phys. 10, 2789 (1980). 

[10] R. Moessner, S. L. Sondhi and P. Chandra, Phys. Rev. Lett. 84 4457 (2000); R. Moessner and S. L. 
Sondhi, Phys. Rev. B 63, 224401 (2001). 



4 



